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Abstract 

REACLIB is one of the most comprehensive and popular astrophysical re- 
action rate libraries. However, its experimentally obtained rates for light 
isotopes still rely mainly on the Caughlan & Fowler (1988) compilation and 
have never been updated despite the progress in many relevant nuclear astro- 
physics experiments. Moreover, due to fitting errors REACLIB is not reliable 
at temperatures lower than 10 7 K . 

In this work we establish the formalism for updating the obsolete 
Caughlan-Fowler experimental rates of REACLIB. Then we use the NACRE 
compilation and results from the LUNA experiments to update some im- 
portant charged-particle induced rates of REACLIB focusing on the proton- 
proton chain. The updated rates (available also in digital form) can now be 
used in the low temperature regime (below \tf K) which was forbidden to 
the old version of REACLIB. 
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I. INTRODUCTION 



Nuclear reactions have long been established as the engines that provide the necessary- 
energy for the stars which use that energy to balance the enormous gravitational pressure 
of the stellar gas(see Ref. [1] and references therein). Every stellar evolution code relies on 
a nuclear reaction rate library which is read by the code before any simulation is performed 
(see for example Ref. [2,6]). The size and accuracy of that library determines the quality 
of the relevant simulation and therefore its nuclear reaction rates should be continually 
improved and upgraded . 

One of the most comprehensive such libraries is REACLIB, updated [5] by the Basel 
nuclear astrophysics group and used extensively in small [4] and large [2] scale simula- 
tions. According to its authors, due to fitting errors, REACLIB cannot be safely used 
at temperatures lower than T 9 = 0.01 (where T 9 measured in 10 9 K) , despite the fact that 
nuclear burning of hydrogen isotopes starts at lower temperatures, such as deuterium burn- 
ing down the Hayashi line. That defect of REACLIB is an undesirable consequence of 
fitting a single formula to an array of data (see below) which spans many orders of mag- 
nitude (sometimes more than thirty!) with respect to a very extended temperature range 
0.01 < T 9 < 10. In fact even nuclear rates which are evaluated close to the upper limit of the 
critical temperature region are expected to be contaminated with similar (although less se- 
vere) fitting inaccuracies. For example, very important studies where the application of the 
old REACLIB rates may be questionable are pre-main and main-sequence stellar evolution 
simulations [4,6] including the solar evolution/neutrino studies where temperatures never 
exceed the value of T 9 = 0.016 (the central value of the present sun is roughly T 9 = 0.0157) . 
Therefore it is now obvious that all stellar evolution simulations which start from a Zero- 
Age Main Sequence star (ZAMS) are forced to apply REACLIB to temperatures well below 
T 9 = 0.01 yielding results which may be inaccurate. 

According to the authors of REACLIB most of its charged-particle reaction rates for 
light nuclei rely on the compilation of Caughlan and Fowler [7] . However, since the publi- 
cation of Ref. [7] there has been a very fertile activity in the field of experimental nuclear 
astrophysics leading to experiments which for the first time reached deep into the most 
effective energy of interaction of astrophysical reactions [8]. New reaction rate compilations 
have appeared either for light nuclei [9] or for heavier ones participating in explosive burn- 
ing [10]. REACLIB has not yet been updated and its light nuclei experimentally-obtained 
rates are obsolete. Moreover, the Caughlan-Fowler [7] rates suffer from another source of 
inaccuracy since many higher energy resonances are lumped into one analytical term which 
is an undesirable oversimplification for plausible reasons. 

The present paper has three objectives: 

a) to use the NACRE compilation (as well as other more recent experimental data) 
in order to partially update the light-isotope experimental reaction rates of REACLIB. 
The update is focused on some of the most important reaction rates of the proton-proton 
chain which (in their updated version) can now be used in the critical temperature region 
T 9 < 0.01. 

b) to improve the REACLIB fitting accuracy in such a way which would allow its 
application to high quality studies related to the destruction of short-lived nuclei in pre- 
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main sequence stellar evolution. 

c) to establish the formalism and techniques which will be used in future more extended 
updates of REACLIB. 

The layout of the paper is as follows: 

In Section II there is a brief description of the formalism used in the evaluation of light- 
nuclei thermonuclear reaction rates. In Section III the main components of REACLIB are 
presented while in Section IV we describe the methodology adopted in order to accomplish 
the above mentioned objectives. In Section V some of the most important reactions of the 
proton-proton chain are updated while the results of the present study are summarized in 
Section VI. 



II. CALCULATION OF THERMONUCLEAR REACTION RATES 

The thermonuclear reaction rate (TTR) for the binary reaction X (a, b) Y is given by 
the formula 

r aX = (l+5 aX y l N a N x {au) (1) 

where N a , N x are the number densities and (au) stands for the reaction rate per pair of 
particles given by the formula: 

E \ 



8 r°° ( E \ 

(au) = J ^ a(E)Eexp(-—)dE (2) 



W (kT) 3/2 

The Kronecker symbol 5 aX takes into account that the interacting nuclei can be identical. 

The cross section a (E) , which appears in the TRR, can be non-resonant or resonant 
according to the range of stellar energies E. If the temperature of the star is such that 
the integrand goes to zero before the cross section strikes a resonance the non resonant 
formalism can be used by adopting the formula 

a(E) = m^e->™ (3) 

where S (E) is a slowly varying function of energy called the astrophysical factor and n is 
the usual Sommerfeld parameter. 

On the other hand if the most effective energy of interaction (see Eq.6) is equal to the 
energy (E r ) of a quasi-stationary state of the ensuing compound nucleus then the cross 
section exhibits resonant behavior which can be described by the Breit-Wigner formula: 

K J k 2 (E — E r ) + T (E) /4 1 ; 

where k is the wave number, (E) and T f (E) are the entrance and exit channel partial 
widths, respectively, T (E) is the total width, and to is the statistical factor given by 
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where Ji, J 2 , J,are the spins of the interacting nuclei and of the resonance, respectively. 

When the Breit-Wigner formula is inserted into Eq. (2) the integrand exhibits maxima 
at E r and at the most effective energy of interaction Eq given by 

E = 0.1220 (Z\Z\A) 1/3 Tg /3 Mev (6) 

where Zi, Z2 are the charge numbers of the interacting nuclei and A is the respective reduced 
mass number A = A1A2 {A\ + A 2 ) 1 . 

The thermonuclear reaction rate (TRR) per pair of particles for an isolated narrow 
resonance E r is given by [1] 



where (wy) Er = u;rjF//r (E r ) 

For light nuclei which capture protons or alpha particles (such as those studied in the 
present paper) the compound nuclei will be produced at low excitation energies where 
the level densities are low. In a such a case the statistical model (i.e. Hauser Feshbach 
model) breaks down [10] and sometimes overestimates the reaction rates by several orders 
of magnitude. Therefore, the total Maxwellian averaged reaction rate N A < au > is 
determined by summing up the contributions of (i) single isolated (narrow) resonances 



N A < au > r = N A g) 3/2 (w7 ) ri T-3/' exp (-§ 



(8) 



and their single non-resonant (tail) contribution 



1/2 A£ „ / ZE f 



N A < au > nr = N A \- j j—^S eff exp (-^ j (9) 
so that 

N A < au> = ^N A < au > n +N A < au > nr (10) 

i 

where we have used the familiar notation for rates used in many popular textbooks and 
articles such as Ref. [1] and Ref. [11] and the index (i) in Eq. (10) indicates a particular 
isolated resonance. Note that S e ff is the effective astrophysical factor which is given as a 
function of the experimentally derived zero-energy astrophysical factor and its derivatives 
(S (0), S' (0)...) by the formula 



'V/ = s (o) 



5 + in(E 0+ ^ kT ) + is "w 



12r S(0) V 36 J 2 S (0) 



El + -E G kT 



(11) 



where the E is the most effective energy of interaction, r is given by r = 4.248 (Z 1 Z 2 AT g ) 1 ^ 3 
and A is reduced mass number: A = A\A 2 {A\ + A 2 ) -1 ). 
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Nuclear astrophysics experiments measure all the components \S (0) , S' (0) ...J of the ef- 
fective astrophysical factor S e ff, the resonance energies E r and the respective partial widths. 
Those data are then inserted into Eq. (10) in order to provide formulas for the thermonu- 
clear reaction rates which are used in stellar evolution simulations. However, thermonuclear 
reaction rate data are more easily used when they are given in tabular forms [5,9] so that 
they can be parameterized into reaction rate libraries by using suitable fitting formulas 
[3,5,10]. Those libraries are then uploaded directly by the simulation code for stellar evo- 
lution and nucleosynthesis calculations. 

III. BRIEF DESCRIPTION OF REACLIB 

REACLIB is a complete library of nuclear reaction rates. Using capital letters 
A, B, C, to denote each isotope (parent and daughter ones) the reaction rate library 
REACLIB would consist of the following components: 

DECAYS 

Beta decays and electron captures: A — ► B 

Photodisintegration and beta delayed neutron emission: A — > B + C 

Inverse triple alpha or beta-delayed two neutron emission: A — > B + C + D 

BINARY REACTIONS 

A + B — >C 

A + B — ► C + D 

A + B — >C+D+E 

A + B — >C + D + E + F 

TRIPLE REACTION 

A + B + C — > D 

A + B + C — > D + E 

Each rate is described by three lines. The first line indicates: a) the participating nuclei, 
b) the source of the reaction, c) the type of reaction (resonant, non-resonant), d) if the rate 
is calculated from the inverse reaction rate or not, e) the Q value of the reaction in MeV. 
The second and third line for each rate give the seven fitting coefficients described below. 

All reaction rates in REACLIB have been derived by using the seven-parametric 
(aj,i = 1...7) fitting formula 

Rtot (ai...a 7 ; T 9 ) = ex V (a x + a 2 T 9 ~ 1 + a 3 T 9 " 1/3 + a 4 T 1/3 + a 5 T 9 + a 6 T 9 5/3 + a 7 lnT 9 )) (12) 

where the reaction rate R tot (ai...af, T 9 ) corresponds to: ln2/t 1 / 2 for decays, N A < 
ab > for binary reactions, N\ < abc > for triple-reactions (N A being Avogadro's 
number), and T 9 is the temperature in units of GK. According to the above men- 
tioned formalism (i.e. Eq. (8) , Eq. (9) , Eq. (10)) REACLIB splits the total charged-particle 
induced rate Rtot into one non-resonant R nr (ct^ r ...a™ r ) and (i) resonant components 
R n (cli ...cli) denoted respectively by the superscripts (nr) and (r). Thus the total R to t 
reaction rate will be 

Rtot = Rnr «...<) + J2 R n (a?...a r /) (13) 
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Note that one doesn't have to include all possible resonances of the compound nucleus 
which is formed during charged-particle capture reactions. It is sufficient to include those 
isolated (narrow) resonances which are relevant to the temperature (energy) regime where 
the rate will be applied to. 

Regarding the use of REACLIB, stellar modelers often apply( [3,4], and references 
therein) REACLIB to the critical region 0.001 < T 9 < 0.01 mentioned in our introduction, 
by using a very limited reaction network of light nuclei. Their decision is partly justified 
by the fact that at such low temperatures there is only a tiny nuclear energy production, 
while as regards nucleosynthesis only very light nuclei are destroyed such as deuterium, 
lithium etc. Sometimes nuclear burning at temperatures Tg < 0.0005 is totally disregarded 
and only decays are considered. It is obvious that the obsolete Caughlan- Fowler [7] rates of 
REACLIB may have been a source of errors to all stellar evolution simulations [2-4] that 
have used it. 

There are various versions of REACLIB which can be downloaded from Ref. [12]. The 
most recent version of REACLIB currently available on-line [12] by the Basel group in- 
volves the unprecedented number of 5.411 isotopes with a mass numbers range 1 < A < 
279. However, its light-isotope charged-particle experimental rates are still those of Ref. [7], 
which underlines the importance of the present study. 

IV. ADOPTED METHODOLOGY 

Fitting the REACLIB fitting formula Eq. (12) to the tabular reaction rate data (e.g. the 
NACRE ones) over the entire temperature range is not the most accurate approach. The 
fitting engine is forced to fit a single formula to an array of data which spans many orders 
of magnitude (sometimes more than thirty!) with respect to a very extended temperature 
range 0.001 < T 9 < 10. This approach generates sometimes a significant error which will be 
avoided in our study. 

On the other hand, there are admittedly more sophisticated mathematical functions 
that could fit the data much better than Eq. (12), such as those given by NACRE [9]. 
However, we must follow the format of REACLIB otherwise we should be dealing with a 
different reaction rate library of a different format (whose fitting function might be taxing 
the computer considerably). 

We must now turn to the format of the NACRE data which must be seriously taken 
into consideration. The NACRE data are given in the form of an array of values (rates 
with respect to temperature plus uncertainties) which is the result of combining various 
individual rate components, namely: non-resonant, narrow resonant+tails, broad resonant 
and multi-resonant rates. It is impossible to extract the individual rates from the combined 
NACRE tabular data although that is necessary in our work due to the format of REACLIB. 
Fortunately, the NACRE authors have also derived analytical approximations to each of 
these rate components, thus providing a tool for uploading the new NACRE rates into 
REACLIB. We use the ORIGIN fitting package which relies on the Levenberg-Marquardt 
(LM) algorithm (one of the most powerful and reliable fitting methods) to perform non- 
linear regression. Actually, we fit Eq. (12) to each of these analytical approximations only 
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over the temperature range where the respective rate component plays a non-negligible role 
to the total rate. Outside this range the new individual REACLIB rates may not be very 
reliable (although the total rate can be safely used). For example, if the non- resonant (NR) 
rate of a particular reaction is considerably smaller than the respective (first resonance) 
Rl rate at temperatures T > T 9 * then our fitting range for the NR rate for that reaction 
would be [0.001, Tg] provided the relevant NR Eq.(12) behaves asymptotically correctly 
at T > Tg (e.g. it is a decreasing function of temperature). In such a case we wouldn't 
recommend using individually the NR rate of that particular reaction at temperatures T > 
T*. However, at T < T 9 * that particular NR rate is suitable for all practical applications, 
where of course the relevant reaction rate obeys the general rule of Eq. (14) 

Normally the fitting procedure would involve fitting Eq. (12) to the NACRE ana- 
lytic functions, however we have decided to fit the exponent of Eq. (12) to their natural 
logarithms , which is a more accurate approach. 

As for the non-resonant rate of resonant reactions, in several cases we have calculated 
the non-resonant rates by using the numerical integration formalism adopted by NACRE 
[9]. This was necessary in order to verify that the analytic formula given by NACRE has 
not been misprinted.. 

Note that in our fit we adopt the assumption made by NACRE [9] that all rates N A < 
uu > which obey the condition 

N A < au >< 10~ 25 (14) 

can be considered negligible in practically all astrophysical applications. 

We assess the quality and relevance of the updated REACLIB by using the following 
three tools: 

Firstly, we can ascertain that the new (updated) total REACLIB rates approximate 
satisfactorily the total NACRE ones by plotting their relative difference RD with respect 
to temperature (Tg): 



TjREACLIB(new) (rp \ -pNACRE (rp \ 

n REACLIB(new) (r p \ 
ttfot {-1-9) 



Secondly, we can assess the relevance of updating the REACLIB rates (i.e. the present 
work) by plotting the variation of the relative difference between the old total REACLIB 
rate and the total NACRE rates, with respect to temperature (Tg): 

jREACLIB(old) (rp X tdNACRE , 



jyn.nA^i.inyoia) (rp \ tdNACRE (rp \ 

, K tot (£9j - tijot ( 9/1 

D RE AC LI B (old) /r p \ 
*Hot \ I 9) 



RD (Tg) = 100^ -MA^ko^ZT, V ~ % (16) 



Thirdly we can assess the deviation between the new and the old REACLIB rates by 
plotting their relative difference with respect to time. Provided that the new REACLIB 
rates approximate well the NACRE ones this tool is also a measure of the relevance of the 
present update: 



t-,RE AC LI B (old) (r p \ D REACLIB(new) (r p \ 
, K tot \ 1 9) — tttot \ L 9) K 

j^RE AC LI B (old) (rr n. 
*Hot l J 9j 



RD (T 9 ) = 100^ ^/cWt^ — % ( 17 ) 
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When necessary, we plot the variation of RD (%) with respect to temperature for two 
different temperature regimes: The first regime is relevant to solar evolution simulations 
while the second one covers the entire temperature range given by NACRE. 

In each figure caption we also include the accuracy (n%) of the analytical approximation 
given by NACRE (see Appendix B of Ref. [9]). Therefore all the updated REACLIB reaction 
rates relying on the NACRE compilation carry an inherent fitting error of (n%) . 

V. UPDATING THE REACTIONS OF THE PROTON-PROTON CHAIN 



This reaction is of paramount importance to stellar evolution (and especially to solar 
evolution/neutrino) studies therefore its thermonuclear reaction rate at relevant tempera- 
tures should be as accurate as possible. 

Adopting the same formula for the astrophysical factor as the one given by NACRE [9] : 



we numerically integrate Eq. (2) and then fit Eq. (12) to the derived tabular data over the 
region 0.001 < T 9 < 0.1. We confined our fit to a much shorter range 0.001 < T 9 < 0.02 but 
the accuracy didn't improve. It is obvious from Fig.l that at temperatures Tg < 0.1 the old 
REACLIB values are very close to the new ones, therefore the old 1 H (p, z/e + ) 2 H rate need 
not be updated as regards the solar evolution zone. However, although at temperatures 
T 9 < 0.1 both the old and the new REACLIB rates approximate well the NACRE data, at 
higher temperatures the old REACLIB rate significantly deviates from the NACRE one. 
This deviation is particularly relevant to explosive hydrogen burning simulations which 
sometimes are performed using the REACLIB library (e.g. the TYCHO [3] code, which 
is based on the REACLIB library, is equipped with explosive burning simulation engines). 
Therefore we recommend the use of the present updated reaction rate parameters over the 
entire spectrum of temperatures: 0.001 < Tg < 10. 

Figure 1. 1 H (p, ve + ) 2 H : The variation RD(%) of the relative difference between the 
REACLIB rates and the values obtained by numerical integration of Eq. (2) using the 
NACRE data. The solid curve represents the RD between the new REACLIB rate and the 
NACRE one while the dotted one represents the RD between the old REACLIB rate and 
the NACRE one. (n = 3%) 



NACRE specifies two temperature regimes and fits two different functions for the re- 
action rates. However, REACLIB cannot follow the same format. Instead, we fitted the 
REACLIB rate formula to the NACRE tabulated rates and we found that the NACRE rates 
and those given by the new REACLIB are almost identical for the entire low-temperature 



A. 1 H(p,ve + ) 2 H 




MeVb 



B. 2 H(p,-/) 3 He 
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regime. However, we observe a 37% maximum relative difference between the old REA- 
CLIB rates and the NACRE ones while at the solar temperature regime, in particular, the 
RD is roughly 25%. All deuterium burning studies which have relied on the old REACLIB 
should take that observation seriously into account. 

Figure 2. 2 H (p, 7) 3 He : The variation of the relative difference RD between the old 
REACLIB rates and the NACRE ones with respect to temperature (n = 3%) 

Figure 3. 2 H (p, 7) 3 He : The variation of the relative difference RD between the new 
REACLIB rates and the NACRE ones with respect to temperature (n = 3%) 

C. 2 H{d,~f) 4 He 

Fitting the REACLIB formula (i.e. Eq. (12)) to the tabulated data of the NACRE 
compilation yields very satisfactory results. The new REACLIB formula fits excellently the 
NACRE tabular data and therefore plotting the variation of the relevant RD with respect to 
temperature is unnecessary. Instead, we plot the variation of the relative difference between 
the old and the new REACLIB rates with respect to temperature. Fig. 4 shows that the 
old REACLIB rate in the solar regime can be up to 15% larger than the rate predicted 
by NACRE while this discrepancy is fixed by the new fit. The updated REACLIB rate 
is practically the same as the NACRE one in the same regime. However, as shown in the 
same figure at larger temperatures the new REACLIB rate (i.e. the NACRE rates) are 
significantly larger than the old ones. 

Figure 4. 2 H (d, 7) 4 He :The variation of the relative difference RD between the new 
and old REACLIB rates with respect to temperature (n = 3%) 

D. 2 H(d,n) 3 He 

The new REACLIB formula fits excellently the NACRE tabular data and thus, as in the 
case of 2 H (d, 7) A He reaction we only plot the variation of the relative difference between 
the old and the new REACLIB rates with respect to temperature. 

According to Fig. 5 in the solar regime the old REACLIB rates are up to 12% smaller 
than the new (or NACRE) ones while at higher temperatures this effect is reversed and the 
old rates become larger than the new ones (up to 90%). 

Figure 5. 2 H (d,n) 3 He :The variation of the relative difference RD between the new 
and old REACLIB rates with respect to temperature (n = 4%) 

E. 2 H(d,p) 3 H 

The new REACLIB formula fits excellently the NACRE tabular data. According to 
Fig. 6 in the solar regime the old REACLIB rates are up to 10% smaller than the new (or 
NACRE) ones while at higher temperatures this effect is reversed and the old rates become 
larger than the new ones (up to 90%) 
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Figure 6. 2 H(d,p) 3 H: The variation of the relative difference RD between the new 
and the old REACLIB rates with respect to temperature, (n = 5%) 

F. 3 He( 3 He,2p) A He 

For this reaction we don't rely on NACRE data in order to produce its REACLIB rate. 
The LUNA collaboration has managed to lower the beam energy of their experiment so 
much that they have recently evaluated the relevant astrophysical factor with the highest 
precision ever. Therefore we follow the most reliable procedure of numerically integrating 
the thermonuclear reaction rate integral, a method followed by NACRE as well. 

Then we fit Eq. (12) to the array of numerical data. The new REACLIB formula 
represents the data very satisfactorily and according to Fig. 7 there is a notable deviation 
from the old REACLIB formula, which may have a non-negligible effect on solar evolution 
simulations using REACLIB. 

In Fig. 7 we plot the variation of the RD between the REACLIB rates (old and new) and 
the rate obtained by numerically integrating Eq. (2) using the most recent LUNA data. 
We observe that the RD between the new REACLIB and the LUNA rates (solid curve) is 
consistently smaller than the respective RD (dotted curve) between the old REACLIB and 
the LUNA rates. Especially in the solar regime the old REACLIB rate deviates from the 
LUNA one by up to 7% whereas, in the same region, the RD between the new REACLIB 
rate and the LUNA one is less than 1%. 

Figure 7. 3 He ( 3 He, 2p) He The variation of the relative difference RD between the 
(old/new) REACLIB rates and that obtained by using the LUNA data with respect to 
temperature. The solid (dotted) curve represents the RD between the new (old) REACLIB 
rate and the LUNA one. 

G. 3 He (a, 7 ) 7 Be 

The new REACLIB formula fits excellently the NACRE tabular data. According to 
Fig. 8 in the solar region the new REACLIB rate approximates the NACRE rate with an 
accuracy of 1% or better, while the old REACLIB rate (see Fig. 9) can be up to 2.5% 
larger than the NACRE one. According to Fig. 10, which shows the deviation between the 
old and the new REACLIB rates, the new REACLIB rate is approximately the same as the 
old REACLIB one. However, due to the importance of the 3 He (a, 7) 7 Be reaction in the 
solar neutrino studies we recommend using the new updated REACLIB rate. 

Figure 8. 3 He (a, 7) 7 Be : The variation of the RD between the new REACLIB rate 
and the NACRE one with respect to temperature, (n = 6%) 

Figure 9. 3 He (a, 7) 7 Be: The variation of the RD between the old REACLIB rate and 
the NACRE one with respect to temperature, (n = 6%) 

Figure 10. 3 He (a, 7) 7 Be :The variation of the RD between the old and the new REA- 
CLIB rates with respect to temperature, {n = 6%) 
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H. 6 Li(p,-f) 7 Be 

The new REACLIB rate approximates much better the NACRE rate than the old one. 
According to Figs. 11 and 12 in the solar region the old REACLIB rate can differ from 
the NACRE rate by up to 80% while the respective discrepancy for the new REACLIB is 
always less than 2%. 

Figure 11. Q Li(p^) 7 Be :The variation of the RD between the old REACLIB rate 
and the NACRE one with respect to temperature, (n = 7%) 

Figure 12. 6 Li (p, 7) 7 Be :The variation of the RD between the new REACLIB rate 
and the NACRE one with respect to temperature, {n = 7%) 

I. 6 Li (p, a) 3 He 

REACLIB distinguishes two components for that rate: a non-resonant and a resonant 
one while NACRE adopts a single non-resonant fit. By fitting Eq. (12) to the single analytic 
formula given by NACRE we observe a very satisfactory representation of all the NACRE 
tabulated data. In Fig. 13 we compare the new REACLIB fit and the old two-component 
one where a minor deviation between the new fit and old one is observed. Accordingly we 
recommend a single non-resonant REACLIB formula for the updated library 

Figure 13. 6 Li (p, a) 3 He :The variation of the RD between the old REACLIB rate and 
the new one with respect to temperature. (2%) 

J. 7 Li (p, 7) s Be 

This reaction is missing from REACLIB and so is the relevant ensuing decay 8 Be — > 
i He + 4 He, therefore we cannot compare the new REACLIB rates to the old ones. The im- 
portance of this reaction to the PPII chain is that it is in competition with the 7 Li (p, a) A He 
reaction, the latter being much more important to the solar evolution studies, of course. 
The s Be produced in the 7 Li (p, 7) s Be reaction, which is unstable and decays into two 
alpha particles in 2.6 x 10~ 16 s, is extremely important to the triple alpha reaction as well. 
Due to the importance of that reaction we will defer its study (and/or update) to a later 
paper where we will also investigate the effects of its absence on the simulations that have 
used REACLIB. 

K. 7 Li (p, a) 4 He 

According to Fig. 14 the non-resonant rate dominates at temperatures Tg < 4. We have 
compared the NACRE rates and those given by the old REACLIB and have found that their 
small differences are within the relevant uncertainties. Therefore, no update was deemed 
necessary for that reaction. 

Figure 14. 7 Li (p, a) 4 He : The logarithms of the NACRE rates (resonant and non- 
resonant) with respect to temperature. (6%) 
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L. 7 Li(a,-f) 11 B 

NACRE evaluates the rates using a non- resonant (NR), a resonant (Rl) and a multi- 
resonant (MR) rate (see Fig. 15) while REACLIB relies only on a NR and a Rl rate. In 
Fig. 16 we plot the variation of RD between the old/new REACLIB rates and the NACRE 
one for all relevant temperatures. We do not include an inset figure for the solar regime 
as that particular reaction is irrelevant to solar evolution studies. It is obvious that the 
updated rates approximate the NACRE ones better than the old REACLIB ones. 

Figure 15. 7 Li (a, 7) 11 B : The logarithms of the NACRE rates (resonant, non- 
resonant and multiresonant) with respect to temperature, (n = 17%) 

Figure 16. 7 Li (a, 7) n B: The variation of the RD between the old (solid 
curve) / new(dotted curve) REACLIB rates and the NACRE one with respect to 
temperature, (n = 17%) 

M. 7 Be(p,-f) 8 B 

The NACRE non-resonant data for this reaction have been superseded by more recent 
ones [13]. According to NACRE [9] the recommended S-factor at zero energy is Sn (0) = 
21 ± 2 eV b while according to Ref. [13] it is should be S 17 (0) = 18.6 ± 1.2 eV b. Despite the 
notable difference in the zero-energy astrophysical factor we decided to use the NACRE 
data for consistency. However, it should be noted that for high quality solar neutrino 
calculations the more recent value [13] should be adopted which would lead to an 11.4% 
decrease in the relevant non-resonant rate. 

The new REACLIB fitting approximates the NACRE tabular data better than the old 
one in the range 0.002 < T 9 < 2.2 while at higher temperatures the old REACLIB rate 
constitutes a better approximation. Due to the large uncertainties involved at such high 
temperatures we recommend the use of the updated REACLIB rates over the entire spec- 
trum of temperatures. In Fig. 17 we observe that the NR component of the rate dominates 
the Rl one everywhere while in Fig. 18 we plot: a) the variation of the RD between the new 
REACLIB rate and the NACRE one (solid curve) and b): the variation of the RD between 
the old and the new REACLIB rates (dotted curve). It is obvious that the new REACLIB 
rates are generally more reliable than the old ones especially at the solar evolution regime. 

Figure 17. 7 Be (p, 7) 8 B :The logarithms of the NACRE rates (resonant, non-resonant 
and multiresonant) with respect to temperature, (n = 3%) 

Figure 18. 7 Be(p,-y) 8 B :The variation of the RD between the new REACLIB rate 
and the NACRE one (solid curve) and the variation of the RD between the old and the 
new REACLIB rates (dotted curve), (n = 3%) 

VI. CONCLUSIONS 

The nuclear reaction rate library REACLIB is one of the most comprehensive and pop- 
ular ones and is extensively used in stellar evolution and nucleosynthesis simulations. In 
the present study some very important light-isotope charged-particle experimental rates 
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of REACLIB have been updated using the NACRE [9] compilation and results from the 
LUNA [8] experiments. We have focused on the most important reaction rates of the 
proton-proton chain while the updated REACLIB rates can be used at temperatures which 
were forbidden to the old ones. The deviation between the new and the old REACLIB 
rates is sometimes significant especially at the low temperature regime (0.001 < T 9 < 0.01) 
of deuterium burning 2 H (d, 7) A He 2 H (d,p) 3 H 2 H (d,n) 3 He (where the old REACLIB 
rates were unreliable). Another notable deviation is that the most important reaction 
1 H (p, ve + ) 2 H appears to be faster in the updated REACLIB than in the old one. The ef- 
fects of these deviations on explosive hydrogen burning and big-bang nucleosynthesis should 
be carefully investigated by adopting successively the old and the new REACLIB libraries 
in relevant simulations (currently under study by the author). 

Another improvement of the new REACLIB rates (which are also available in the same 
digital form like the old ones) is that we have improved their fitting accuracy in such 
a way which allows their application to pre-main sequence stellar evolution. Finally we 
have established the formalism and techniques which will be used in future more extended 
updates of REACLIB (soon to appear by the author). 
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REACTION 


NR/R 


al 


a2 


a3 


a4 


a5 


a6 


a7 


'H(p,ve + ) 2 H 


NR 


U . 3894ZUE+UZ 


U . UU 1 b bUE + U U 


U . UU 1 b bUE + U U 


U . 89 / 3 4 UE+ U 1 


U.04U1/UE+UU 


U . Z 983U UE— Ul 


u . z z id/ yE+u i 


2 H{p,y) 3 He 


NR 


. 635671E + 01 


-0.001600E-01 


-0.368571E+01 


. 328570E+01 


-0.201860E+00 


. 129900E-01 


-0.761200E+00 


2 H(d,y) 4 He 


NR 


0.264093E+01 


0.053000E-02 


-0.448526E+01 


. 146128E+01 


. 549700E+00 


-0 . 697900E-01 


-0 . 105681E+01 


2 H(d,n) 3 He 


NR 


. 184878E + 02 


0.002200E-01 


-0.438956E+01 


0.246392E+01 


-0 . lb6590E + 00 


. 097600E-01 


-0 .099180E+01 


2 H(d,p) 3 H 


NR 


. 189686E+02 


0.017000E-02 


-0.435280E+01 


. lb5054E+01 


-0 . 012120E + 00 


-0 . 519000E-02 


-0 .890840E+00 


3 He( 3 He,2p) 4 He 


NR 


0.287852E+02 


-0 .238000E-02 


-0.114354E+02 


-0 . 59109bE + 01 


. 103574E+01 


-0 . 076820E + 00 


. 6999bOE+00 


3 He ( 4 He, y) 'Be 


NR 


U . 1 / Z U / UE + UZ 


-U . UU lb / UE + U U 


— U . 1Z443ZE+UZ 


— U.Zob/9/E + Ul 


U . Z4138UE+UU 


p p n o^pptp^pp 
-U . U13oUUE+UU 


n p r p p i ptp ^ p p 
-U . Ub9UiUE + UU 


6 Li{p,yy Be 


NR 


. 156649E+02 


-0.420000E-03 


-0 . 820981E+01 


-0 . 240965E+01 


. 324690E+00 


-0 . 021530E + 00 


-0 .230010E+00 


6 Li(p,a) 3 He 


NR 


. 253677E + 02 


-0.430000E-03 


-0 . 823640E+01 


-0 . 152610E+01 


. 144720E+00 


-0 . 007170E + 00 


-0 .331910E+00 


7 Li{p,y) "Be 


NR 


* 


* 




* 


* 


* 




7 Li(p,a) 4 He 


NR/R 


* 






* 


* 


* 




7 Li(a,y) u B 


NR 


. 642203E+03 


-0 . 143030E+00 


. 545753E+02 


-0 . 9b5712E+03 


. 4b0363E+03 


-0 . 222558E+03 


. 1617b9E+03 


7 Li(a,y) u B 


Rl 


. 581419E+01 


-0.295900E+01 


0.403880E-08 


-0 . 516330E-08 


0.2b5870E-09 


-0 . 140696E-10 


-0 . 150002E+01 


7 Li(a,y) u B 


R2 


. 924966E+01 


-0 . 492205E+00 


-0.9303bOE-08 


. 122013E-07 


-0 . b491bE-09 


. 354410E-10 


-0.023000E+00 


7 Be{p,y)«B 


NR 


. 124677E+02 


0.370150E-06 


-0 . 102b43E+02 


. 863000E-02 


-0 . 06523E+00 


0.046720E+00 


-0 . 667580E+00 


7 Be{p,y)«B 


Rl 


. 762567E+01 


-0.734500E+01 


-0 . 10bb70E-07 


. 131810E-07 


-0 . bb3960E-09 


. 347160E-10 


-0 . 150002E+01 



Table 1. The updated values of the seven fitting parameters a,...a 7 used by the REACLIB reaction rate library. The asterisks indicate values 

that have not been updated for reasons given in the text while the notation NR/R indicates: NR=Non-Resonant reate, R=Resonant rate. A 
digital file with the updated rates is available which can be readily incorporated into the REACLIB library (available upon request). 



